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ABSTRACT 

Using analytical models and cosmological N-body simulations, we study the free- 
free radio emission from ionized gas in clusters and groups of galaxies. The results 
obtained with the simulations are compared with analytical predictions based on the 
mass function and scaling relations. Earlier works based on analytical models have 
shown that the average free-free signal from small haloes (galaxies) during and after 
the reionization time could be detected with future experiments as a distortion of 
the CMB spectrum at low frequencies {v < b GHz). We focus on the period after 
the reionization time (from redshift z = up to z = 7) and on haloes that are more 
massive than in previous works (groups and clusters) . We show how the average signal 
from haloes with M > IQ^^h^^MQ is less than 10% the signal from the more abundant 
and colder smaller mass haloes. However, the individual signal from the massive haloes 
could be detected with future experiments opening the door for a new window to study 
the intracluster medium. 



1 INTRODUCTION 



In the near future, new high sensitivity experiments 
observing at radio and millimeter wavelengths will open a 
new window to study the high redshift Universe and in par- 
ticular the re-ionization period. Among these exper iments, 
the Square Kilometer Array (or SKA hereafter, iTavloJ 
|2000.) and the Atacama Large Millimeter Array (or ALMA 
hereafter, iKurz fc Shaveilll99dl ') are the most relevant ones 
due to their sensitivity and angular resolution. These 
experiments will be able, for the first time, to trace in detail 
the distribution of neutral hydro gen before re-ioniz ation 
(through the 21 cm line, see e.g. iSchneider et al.ll20oj 'l. and 
the transition between a neutral and ionized Universe at 
the time of reionization (ALMA could see the first galaxies 
emerging at the reionization time). 

The study of the re-ionization period will offer a unique 
window to help us understand the formation of the first stars 
and galaxies. The possibilities of this new window for astron- 
omy has motivated many studies that focus, for instance, 
on th e 21-cm line radiation from neut r al gas ijScott fc ReesI 
1990'; 'Miralda-Escude & Rees 'l99^; 'Kunth ot a ll 1 19981: 
Oh fc Ma ck 2003; Tozzi et al., 2000; Ciardi fc Madaul 

20041 : 



Furlanetto et 'al. | 20041: Zaldarriaga et al 



Burigana et al.ll2004l : lGnedin fc Pradall2004l'), or on the kine- 
matic Sunyaev-Zeldovich effects (kSZ. ISunvaev fc Zeldovi"chl 
198d ') from inhomogeneous (patchy) reionization on large 



scale s (|Santos et all l2003l : llliev et al.l l2007l : Ijelic et all 
l2010t ). In lOhl (| 19991 '). it is proposed that the reionization 
can be studied also through the H-a emission, useful to 
trace young star formation regions. 

Another signal emerging from the ionized regions will 
be the free-free from interactions between the electrons and 
ions in the plasma. The photons emerging from these in- 
teractions can be observed in the radio and microwave 
bands. The distortion that free-free induces on the back- 
ground temperature in the Rayleigh- Jeans part of the spec- 
trum (Bartlctt fc Stcbbins 1991) is ac tually co nstrained by 
the ground based measurement of lBersanelli et al.l (|l994l ) at 
2 GHz, Uff < 1.9 X 10-"^ (95% CL). 

Most of the efforts focus on the study of the 21-cm 
line and the interaction between the CMB photons and 
the ionized clouds but little has been done in relation to 
the free-free signal. In this paper we focus on the free-free 
emission and its ability to trace the ionized medium. The 
free-free emission (or bremsstrahlung) can be potentially 
observed in the local Universe and up to the re-ionization 
era. UV radiation emerging from the first stars and quasars 
ionized the neutral hydrogen creating expanding bubbles 
of ionized plasma. During a free-free interaction between 
two charged particles (free electrons and ionized atoms), 
the electron loses part of its kinetic energy by emitting 
a photon. The energy of the photon ranges from the 
radio to the X-ray wavelength depending on the electron 
temperature. Since this interaction involves two particles. 
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its intensity depends on the square of the free electron 
(or equivalently the ion) density, ric- This dependence 
makes the free-free signal an interesting candidate for 
cross-correlations with other signals like the SZ effect where 
the signal amplitude depends linearly on rie. 

In the late 90's, an experiment was designed to measure 
the distortion of the CMB spectrum due to the cosmological 
free-free signal; the Absolute Radiometer for Cosmology, 
Astrophysics and Diffuse Em i ssion (or ARCADE, see 
Fixsen et aLlj2004 iKogut et all l2006l : iFixsen et al] l2009l : 
Seiffert et al.ll201ll for details). Its goal is to detect the aver- 
age free- free signal at frequencies around 1 GHz. Studying 
the distortion of the CMB spectrum at these frequencies 
would allow, in principle, to set strong constraints on 
the history of reionization of the Universe. Recently, the 
ARCADE team presented the results of the ARCADE2 
mission that studies bot h Galactic and extragalactic 
signals ([Fixsen et al.l [200^ . They detect a signal that is 
significantly larger than the expected extragalactic radio 
background (a factor ~ 5 brighter than the estimated 
contribution from radio point sources). The ARCADE 
team is currently exploring the possible causes of such a 
signal like for instance possible foreground contamination, 
synchrotron emission from Earth's magnetic field or CII 
lines. In the latest review of the results of the mission 
l|Seiffert et al.|[201ll ). the authors still report that the excess 
detection remains unexplained, even though the three 
main sources of errors. Galactic emission, instrumental 
systematic errors and radio emission from the faint end of 
the distr ibution o f know n sources, are carefully taken into 
account. ISharpel (|2009l ) has suggested that the observed 
excess is produced by synchrotron radiation emerging from 
an optically thin low density magnetized plasma region in 
the heliosphere of the Sun. 

[OhI l|l999f) presents an exhaustive treatment of the dif- 
ferent sources of radiation that could be detected with SKA 
and ALMA in the range of the radio frequencies. He pays 
special attention to the free-free signal from small haloes 
and the intergalactic medium (or IGM) and concludes that 
the IGM signal is subdominant when compared with the 
sig nal from haloes. Ano t her in teresting work is presented 
in ICoorav fc Furlanettd (|2004 ) where the authors use a 
halo model to predict the amount of free-free signal from 
haloes. The authors also compute the angular power spec- 
trum of the signal produced by the free-free below 2 GHz 
and make prediction s in th e context of the ARCADE mis- 
sion. IBurigana et al.l l| 19951 ) discusses different physical pro- 
cesses involving the CMB photons a nd the ionized medium , 
including also the free-free signal. In IBurigana et al.l (|2004l ') 
the authors discuss about the possibility of detecting the 
individual sources of free-free signal. 

All these works have focused on the signal from small 
and cold haloes, largely ignoring the signal coming from 
larger and hotter haloes (groups and clusters) . In this paper 
we will study the regime of more massive haloes and focus on 
the period after reionization. Also, an important advantage 
of working with more massive haloes is that their modeling 
is much simpler than in the case of smaller haloes. The cool- 
ing time is significantly larger for massive haloes and one can 



more easily ignore highly non-linear phenomena like radia- 
tive cooling. 



2 FREE-FREE EMISSION 

In a hot plasma with temperature T, the electrons 
move with kinetic energy Ee = 3/2khT where kh is 
the Boltzmann's constant (1.38 x 10"^^ J/K). The min- 
imum T required to i onize a plasma is ~ 2 x 10* K 
ijZaldarriaga et al]|200ll ). This is also the temperature at 
which most of the cooling rad iation occurs in a typical 
galaxy halo (|Fardal et al.|[200ll ). This temperature can be 
translated into a kinetic energy for the free electrons, typ- 
ically in the order of 2 x 10^^ erg (~ 1 eV). The collisions 
between opposite charged particles within the plasma mod- 
ify the paths of the electrons that lose a few percent of their 
kinetic energy {bremsstrahlung or brake radiation). The net 
effect is the bulk emission of photons in the radio frequency 
range (1 ~ 10 GHz), strongly dependent on the square of 
the electron density. Note that this square dependence im- 
plies a crucial role of the density contrast pattern inside the 
haloes. 

The bremsstrahlung, or free-free signal, can be 
parametrized in terms of the el ectron density rie and tem- 
perature Te as (see for example 'Rvbicki fc Lightmanlll979l : 
ICoorav fc FurTanetto 2004 ; Oh fc Mack. 2061 ): 

6, = 5.4 X IQ-^^n^T-i/^ffff (:., Te) exp , (1) 

in units of ergs cm~'^ s~^ Hz~^ The Gaunt factor 

(iKarzas fc Latter] Il96ll : IBurigana etHI Il995l ). Qff, IS com- 
puted for the observed frequency but it has a weak depen- 
dency on the temperature of the gas. In the Rayleigh- Jeans 
limit (where the free-free radiation is more relevant) the ex- 
ponential part can be safely neglected. From Eq. U]) it is 
clear that the free-free emissivity depends mostly on the 
electron density no. The inverse dependence with the root 
square of temperature is a direct consequence of the thermal 
Maxwellian distribution. 

The luminosity of an ionized volume of space with con- 
stant Tie and Te can be obtained from Eq. ([1} by inte- 
grating the electron density and temperature over that vol- 
ume. By dividing this luminosity by the corresponding lumi- 
nosity distance, the flux (or brightness) in Jy can be derived 
(1 Jy — 10~^^ergs/s cm'^ Hz). In terms of the temperature 
distortion, the brightness can be transformed into equivalent 
temperature. In the Rayleigh- Jeans limit we have that 

AT oc FX^ (2) 

where F is the free-free flux and A — cjv. Thus, at 1 GHz, 
while the flux does not change much with frequency, the 
free-free temperature distortion is expected to be higher 
than at 10 GHz by a factor 100. This simple scaling shows 
the convenience of looking for the free-free signal at lower 
frequencies. Several attempts have been made in the past 
to measure the free-free distortion at low frequencies as 
a deviation of the nearly perfect CMB blackbody energy 
spectrum. The first accurate measurements of the spectrum 
of the CMB were made by the FIRAS and have shown 
no departure from the blackbody spectrum (within the 
error bars) in the frequency range of 60 — 600 GHz. It 
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is expected that new experiments will detect the average 
free-free contribution at lower frequencies in the shape of a 
distortion of the CMB energy spectrum. 



3 PREDICTIONS FROM ANALYTICAL 
MODELS 

Through analytica l hal o mod els it is possible to explore 
a wide range of cases. lOhl j 19991 ) shows that the free-free 
contribution coming from the diffuse IGM is significantly 
smaller than the signal from ionized haloes so it can be safely 
ignored. Two ingredients are needed in order to compute the 
average free-free signal from haloes. First, the mass function, 
n{M, z), that predicts the average number of haloes per red- 
shift, z, and mass, M, intervals and, second, a model for 
the internal gas distribution (and temperature) inside the 
haloes. The abundance of haloes can be computed from the 
mass functi on given a cosmol o gical model. We use the mass 
function of ISheth fc TormenI (| 19991 ) (or ST mass function 
hereafter) for this purpose. The ST mass function repro- 
duces well the results obtained with large N-body simula- 
tions. For the internal distribution of the gas in the haloes 
and temperature we assume a standard isothermal /3-model 
with /3 — 2/3. The gas density profile plays an important 
role since steeper profiles can produce a larger free-free sig- 
nal with the same amount of gas (as it happens in the X- 
ray band with gas in galaxy clusters). Other more realis- 
tic models can be found in the literature (see for instance 
lAscasibar et al.ll2003l : lAscasibar fc Diegoll2008l ) but for sim- 
plicity we will use the /3-model as this model requires only 
three parameters (To for the temperature and -Rc, and Uo 
for the /3-model). 



3.1 Predictions for a single halo 

The /3-model is widely used in the context of 
galaxy clusters to des cribe the electron density as a 
function of radius (e.g. ICavaliere fc Fusco-Femiaiicl 1 19761 : 
iDiego fc Maiumdarll200i ): 



■ R 



(3) 



where no is the electron density at the center of the halo 
and Rc is the core radius and we have assumed = 2/3. 

For the /J-model, the free-free luminosity can be com- 
puted integrating Eq. ([1} over the volume of a sphere of 
radius i?vir 



e,dV 



= Cnlni ^tan ^ ^^^^^ 



(4) 



where C = 5.4 x 10"^^27rTe ^^'^gg{v,Tc) and the argument 
p is the ratio Rvii/Rc- For simplicity we have dropped the 
negligible term exp{—hv / {khTe)) in Eq. ([T]). 

The halo luminosity can be transformed into flux given 
the luminosity distance, Dl, from the halo at redshift z to 
the observer (at 2 = 0). 



ATvDLiz)'^ 



(5) 



The values of rio, -Rc and T can be computed from scal- 
ing relations. In order to establish these relationships we 
assume that the total mass of the halo, M, and the total 
mass of the gas, Mgas, are proportional to each other with 
the proportionality constant being the universal baryon frac- 
tion, /b = Mgas/M. Given a viria l mass for the halo, the 
virial radius can be expressed as (jFinoeuenov et all 1200 ll : 
IVerde et al]|200ll ) 



i?vir = 1.3Mi/^(l +2 



(6) 



where M15 is the halo total mass expressed in 10^^ h^^ Mq 
and the radius is scaled with the expansion factor (1-1- 
z)~^. Within the virial radius, the relation between the bary- 
onic mass and the electron density profile (given in Eq. ([3])) 



Mg 



fimpnc{R)dV 



(7) 



47r/imp 



-Rvir 



nc{R')R'^dR'. 



Then, a relationship between rio and the total mass of 
virialized halo, Mv ~ Mgas//b, can be established: 



Mv/b 



{p — tan 1 {p))4-Kfj,mpR^ ' 



(8) 



where p — Rvh/Rc- The ratio between _Rvir and _Rc is as- 
sumed to be constant (Rvii/Rc = 10). For /b we assume 
/b ~ 0.13. Finally, for the tem perature we use the relation 
obtained bv lDiego et al.l l|200lf ) which was shown to be con- 
sistent with X-ray measurements 



T(keV) = lOMf/^ 



(9) 



Once no, -Rv, Rc and T are known, it is possible to 
compute the total free-free luminosity, flux and temperature 
distortion of the halo at redshift z from equations Q, ([Sjl 
and ([2]) respectively. 

3.2 The abundance of haloes: the mass function 

In this wo rk we use Sheth fc To rmen (ST hereafter) 
mass function (jSheth fc Tormenll 19991 . see Appendix A). 

We compute the mass function between the masses 
10** ^ M ^ lO^'^h'^MQ. This mass interval covers the 
range from the smallest ionized haloes to the largest galaxy 
clusters. We include the small haloes in our calculation 
for comparison purposes with earlier works and with the 
more massive haloes. In our simple model we will make the 
assumption that all haloes included in the mass function are 
fully ionized. This is not properly true in the low end of the 
mass interval since, as it was discussed in Oh- (199^), the 
low mass haloes will stay ionized only for a limited amount 
of time. Consequently, at a given redshift, only a fraction 
of the low mass haloes are active or fully ionized. The 
conclusions derived from our calculations should be then 
considered as an upper limit in the low mass range (M 
^ 10^^ Mq). On the other hand, the high mass end 
haloes can be considered as fully ionized as most the gas in 
these massive haloes (clusters) will remain ionized by the 
high temperatures of the plasma in the clusters. Regarding 
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Figure 1. Mass function for different redsfiift intervals. Note flow 
small haloes are common at all redshifts (their population drops 
at redshifts larger than the ones shown here). 

the redshift range we will consider only the contributions 
up to re d shift z = 7. The r eionization peri od wa s studied 
in iLoebl (|l996l ): lohl (|l999f ): [oh fc Macd l|2003h . In Fig. 
[T] we show how the mass function behaves for different 
redshifts. The low mass haloes (~ 10* solar masses) show a 
nearly constant abundance at all redshifts while the number 
of massive haloes decreases with redshift. 

In the next subsection we will combine the predicted 
flux of the /3-model from one halo with the mass function 
to compute the mean free-free signal from a cosmological 
volume. 



3.3 Average free-free emission from an analytical 
halo model 

Combining a model for the gas distribution inside a 
halo, like the /3-model, with the abundance of haloes as a 
function of mass and redshift, it is possible to compute the 
mean free- free signal in a solid angle as a function of redshift 
and/or mass. We can also integrate this information in the 
redshift-mass space and estimate the mean free-free signal 
from all these haloes. 

Given a redshift and mass interval, we compute the 
number of haloes in the interval and compute the free-free 
flux from those haloes. After integrating over the entire red- 
shift range (0 < 2 < 7) and mass range (10* < M < 10^^') 
we compute the mean free-free flux from all the haloes. The 
flux is converted into thermodynamic temperature to 
compute the AT/T as a function of the frequency. The 
resulting distortion from our analytical model is shown in 
Fig . |2] When comparing our results with those obtained by 
lOhl ()l999l ). we find th at our mo del (solid line) falls below 
the predicted value by|o3 (| 19991 ). This can be explained by 
the fact that we are assuming a higher temperature for the 
haloe s. Fixing the temperature to T = 10'* K (l i ke in lOhl 
ligggi ). our model predicts mor e signal th an in [Ohi j 19991 ). A 
possible explanation is that, in [Ohl (|l999l ). only a fraction of 
the haloes were active while in our case all haloes are ionized. 

It is interesting to show how the free- free signal depends 



Figure 2. Average temperature distortion due to free-free as 
a function of frequency for our analytical model. The solid line 
shows the distortion obtained assuming that the temperature of 
each halo was computed with the scaling law r(keV) = lOM^^^; 
for reference, we show as a dotted line the distortion correspond- 
ing to a fixed temperature of lO** K for all haloes. The dot-dashed 
line represents the 95% con fidence level observational upper limit 
derived from l lBersanelli et al. 1994). The star and the triangle 
represent the results from (Oh 1999) related to the diffuse IGM 
(AT = 6.0 X 10-'^ K) and to point sources (AT = 3.4 x lO'^ K) 
respectively. 

on the redshift and the mass of the haloes. In Fig. |3]we show 
the free-free signal for different mass intervals. In each inter- 
val, we compute the mean free-free distortion as a function 
of redshift. The smaller but more abundant haloes give a 
larger signal. Also, as we show earlier, smaller haloes have 
more or less the same abundance at all redshifts and hence 
their average free-free contribution shows a slow dependence 
with redshift. Note how the simulation predicts significantly 
less average signal than the analytical model. As we will see 
later, this is a direct consequence of the lack of resolution in 
the simulation that is not able to capture the contribution 
from the smallest haloes. 

This prediction, however, should be taken with care 
since on one hand we assumed that haloes remain ionized at 
all times and the temperature of the gas corresponds to the 
virial temperature of the halo. In small systems, the cool- 
ing time is short and the gas can cool down significantly, 
become neutral and form stars. Our assumptions are only 
valid for the most massive haloes (groups and clusters) and 
the model predictions are robust only in that regime. For 
these objects, the average free-free distortion is of the order 
of a few to several tens of /iK at 1 GHz. Also, the plot shows 
the average free-free signal obtained from simulations (see 
below) as a function of redshift (stars). As we will see later, 
the smaller range of halo masses of the simulation predicts 
a smaller average free-free signal. 

In Fig. |4] we show more explicitly the dependency 
of the average free-free distortion with the mass range 
but for different redshift intervals. Again, smaller haloes 
contribute more to the average signal than massive ones at 
all redshifts. This result is streng thened by the the model 
proposed in iMiniati et al. I (|2004l ). The authors, referring 
to the component of the UV luminosity produced by the 
thermal emission from gas accreting on to dark matter 
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Figure 3. Free-free signal for different mass intervals as a func- 
tion of redshift and for u = 1 GHz. The points represent the 
distortion AT computed from a cosmological simulation of 300 
Megaparsecs (see Sect. |4] below). The minimum mass resolved in 
the simulation at z = 1.57 is M^i^ .^^^ = 1.14 X lO-*^-*^ and the 
maximum mass is iV/max,sim — 

8.21 X 10". 



haloes, calculate that the most contribution is produced by 
haloes with temperatures between 10^ K and a few x lO'^ 
K, corresponding to masses 10" - 10" solar masses. 

It is important to note, though, that cooling and star 
formation play a critical role in determining the actual con- 
tribution of galaxy-sized haloes (M < 10^^ Mq) to the 
temperature distortion of the CMB. On the one hand, the 
temperature of the ionized gas will be around 10* K, much 
lower than predicted by Eq. (|9]), and its density will be con- 
siderably higher than predicted by the /3-model. The com- 
bined action of both effects can boost the expected free-free 
signal by a large factor. On the other hand, a significant frac- 
tion of the gas will be transformed into stars and most of 
the interstellar medium will be in neutral rather than ionized 
form, and therefore it will not emit any bremsstr ahlung ra - 
diation. The net effect is diflicult to quantify, and[o3 (|l999l l 
has resorted to a phenomenological parameter describing the 
fraction of active galaxies or, equivalently, the average ion- 
ization fraction of the gas. These works focus on the signal 
from small and cold systems, where UV radiation from stars 
and quasars ionizes the surrounding neutral hydrogen and 
creates expanding bubbles of ionized plasma. 

Gas cooling, star formation, and feedback processes de- 
termine the amount of ionized gas, its characteristic density, 
its temperature, and thus the total bremsstrahlung luminos- 
ity. We discuss in more detail the regime of hotter, more 
massive objects, where the gas is heated collisionally rather 
than photo-ionized. These systems contribute only to a small 
fraction of the overall cosmological signal. In this work, we 
provide robust lower limits for the signal produced by mas- 
sive objects, based on a simply physical modeling. A more 
detailed treatment of cooling and photo-ionization of the 
interstellar medium will be addressed in a future work. 
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Figure 4. Free-free signal as a function of the halo mass and for 
different redshift intervals. 

4 N-BODY SIMULATIONS 

In the previous section we have shown how the average 
free-free signal from haloes depends on their redshift and 
mass distributions. We also discussed how these predictions 
depend on the internal gas distribution inside the haloes. In 
this section we use numerical simulations to compute the 
free-free signal. Through N-body simulations we can obtain 
the distribution of the electron density, its temperature and 
ultimately the free-free effect which can be projected into 
sky maps. 

We use the GADGET-2 code llSpringellboOSl ). The code 
is a combination of a Particle Mesh Refinement algo- 
rith m and the Tre eSPH method by Hernquist and Katz 
(Hcr nguist fc Kat3[l989i ). For the cosmological parameters 
we use the concordance model: Qa ~ 0.73, = 0.27, Qh ~ 
0.039, = 0, (78 = 0.79, h = Ho/{WO Km s-^Mpc-i) = 
0.72 where erg is the RMS mass fluctuation on a sphere of a 
radius of 8 Mpc. 

We create th e initial conditions at redshift z=49 with 
the code 2LPT (|Crocce et al.l |2006| .) based on a second- 
order Lagrangian perturbation theory. The initial condition 
is evolved with GADGET-2 from z; = 49 until z = 0. For the 
main simulation, we use a cosmological volume with 512'' 
particles of dark matter and 512^ particles of gas distributed 
in a box size of (300/i~^ Mpc The force smoothing pa- 
rameters has been set to 1/30 of the inter particle distance, 
and corresponds to 20 Kpc for the 300/i~^ Mpc simulation. 

The outputs (or snapshots using the GADGET-2 termi- 
nology) of the 300/!.~^Mpc box were chosen at redshifts for 
which the comoving distance between both ends of the box 
would overlap between consecutive redshifts. Each snap- 
shot is analyzed independently from the others. We as- 
sume that the Universe is fully ionized below z=7 and we 
concentrate on this regime. The masses for the dark mat- 
ter and gas particles are Mdm = 1.165 X 10^°/i"^Mq and 
Afgas = 0.17 X 10^°/i"^Mq respectively. 

The minimum and maximum masses of the structures 
found in our simulation depend on the simulated volume, 
particle mass, and of course the redshift. As discussed ear- 
lier, the free-free signal has a wide dynamical range in mass. 

The choice of the comoving volume of the simulation 
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box is important: on one hand we want to have the largest 
possible box so we can include more massive haloes, on the 
other hand, the small structures have a very significant im- 
pact on the average free-free signal and is also important to 
capture the small scale signal. 

To address the issue of resolution in the N-body simu- 
lation we make a different simulation (with the same cos- 
mology) but with higher resolution. The use of different 
box sizes and resolutions is useful to study a wider range 
of mass es (or resolutions) with N-body simulations (see fo r 
instance lRefregier fc Tevssiej|2002l : iTrenti fc Stiavellill2008l ) . 

The box size of the second simulation is (50/i~^ Mpc)'^ 
(that is 6'^ times smaller in volume). We maintain the same 
number of particles (512^ for dark matter and 512^ for 
gcis). The resulting dark matter and gas particle masses are 
Mdm = 6.1 X IQ'h-^MQ and Mgas = 0.873 x lO'^fe-^M© 
respectively. The masses of the particles are proportional to 
the volume of the simulation boxes divided by the number 
of particles, that is, since the number of particles is the 
same in both simulations, the particle masses are 6'' times 
larger in the 300 Mpc box than in the 50 Mpc one. 

Our simulations do not include cooling nor radiative 
transfer. In a future work we plan to include these mech- 
anisms and improve the predictions. We also plan to ex- 
tend the redshift range into the reionization period. For the 
present work, our intention is to explore the redshift range 
< 2 < 7 and focus on the most massive haloes for which 
the above effects are not so relevant. 



4.1 Range of halo masses in the N-body 
simulation 

We use a halo finder to map the distribution of haloes 
in mass and to associate each simulations with a proper free- 
free emissivity mass range (Fig. [3] and Fig. U)). 

In order to identify haloes and subhaloes in our simu- 
lations we have run the MPl-|-OpenMP hybrid halo finder 
AHF0. A detailed description of AHF is give n in the code 
description paper (jKnoUmann fc Knebell2009l ). We provide 
a brief summary of the mode of operation. By virtue of the 
adaptive mesh hierarchy employed to sample the density 
field, AHF locates overdensities as prospective halo cen- 
ters. The local potential minima are computed for each 
of these density peaks and, treating the prospective halo 
in isolation, the gravitationally bound particles are deter- 
mined. Only peaks with at least 20 bound particles are con- 
sidered as haloes and retained for further analysis. For each 
halo, we compute the virial radius rvir, that is the radius 
r at which the density M(< r) — (47rr^/3) drops below 
Avirp. Here p is the cosmological background density. The 
threshold Avir is computed using the spherical top-hat col- 
lapse model and is a function of both cosmological model 
and time. 

Applying the AHF to the 300 Mpc simulation, we have 
found that the mass of the inside haloes only span between 
M « 10" /i"^M0 and M ^ lO"/i"^M0. Low mass haloes 
are not present in the simulation due to the resolution. On 



^ AMIGA halo finder, to be downloaded freely from 
|http : //wwM.popla. f t .uam. es/AMIGA| 



the high end mass, the limited volume of our simulation 
prevents us from having the most massive clusters in our 
simulation. 



4.2 Average free-free from the simulation 

For each gas particle in the volume, we compute the 
free-free luminosity assuming that the electron density is ap- 
proximately constant over the volume of the particle. Then, 
the integral of the squared of the electron density can be 
computed as: 



rLcdV — Tic — - — . 



finip 



(10) 



iv Jv 

The gas density, Wo, at the position of the particle is 
extracted from GADGET-2 and then transformed into conve- 
nient cm~'^ units. 

Eq. (|10p is used to compute the particle luminosity from 
Eq. U]). The particle luminosity is transformed into particle 
flux given the luminosity distance, Dl, from the particle at 
redshift z to the observer (at z — 0). 



SsiJy) 



(11) 



The internal energy is given by GADGET-2 in units of 
fkm sec~^l^ which is converted into K with the factor: 



CK = 10''(7-l)Af 



(12) 



where 7 = 5/3 is the adiabatic index for a monoatomic ideal 
gas. The scale factor 10*' accounts for GADGET-2 's internal 
units, and kh is the Boltzmann constant. 

After the flux per particle is computed, the fluxes 
are projected along the line of sight into a pixelized 2D 
map. Since the apparent angular size of each box depends on 
the redshift, we have to restrict our analysis to the smallest 
field of view that in our case corresponds to the apparent 
size of the most distant box (about 3 degrees for the 300 

Mpc box). Because we want to compare the distortion 
that our model induces on the CMB temperature as a func- 
tion of frequency, we extract a mean flux, Sa, from all the 
projected maps. The resulting mean brightness is converted 
into temperature in K (antenna or thermodynamic, since we 
are considering low frequencies): 



ATiu) 



cSff 



(13) 



In Fig. [5] we show the result obtained from Eq. (|13p . 
We plot, based on the actual constraint of 1% {dashed line), 
the corresponding upper limit distortion produced by the 
free-free emission over the CMB temperature. The solid line 
shows the mean temperature distortion of the projected map 
in the sky from the simulation while the dot-dashed line 
refers to the assumption of a const ant temperature for all 
particles of T = 10* K fsee [Ohlll999l '). Because in Eq. P the 
strength of the signal depends inversely on the temperature, 
the lower the temperature of the gas, the higher the signal. 

From Fig. [S] we conclude that the average free-free dis- 
tortion predicted from our 300 h"^ Mpc N-Body is weU 
below the current observational constraint {dashed line). A 
much lower temperature for the gas (10* K) in Eq. ([T]) does 
not change the effect too much showing the anticipated weak 
dependency of the free-free distortion with the temperature. 
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Mpc simulation, the smaller box produces a larger free- 
free signal. Considering a slice of 50 Mpc at redshift 
1.57 in both cases; in the 300 Mpc box the average AT 
is AT ^ 10~® K at 1 GHz while in the 50 h'^Mpc box 
AT ^ 5 X 10"*^ K also at 1 GHz. This is a factor 5 more 
signal in the higher resolution case. As shown earlier, this 
extra signal comes from the lower mass haloes although we 
should keep in mind that in our model we are not including 
neither radiative cooling nor partial ionization of the low 
mass haloes. These effects compensate each other partially 
(in terms of the free-free signal) but they will change the 
amount of free-free predicted by our model (again, in the 
low mass haloes more than in the massive ones). 



Figure 5. Free-free emission distortion from the 300 /i-l Mpc 
simulation. The emission is computed from 20 snapshots within 
the interval < z < 7. The field of view covers 2.7 degrees. The 
solid hne shows the distortion in the case where the emissivity 
has been computed with the temperature derived from the sim- 
ulation. The dash-dotted line corresponds to the case where the 
temperature for all particles has been fixed to 10"* K. The dashed 
line s hows the observationa l upper limit constraints (95% CL) 
from l lBersanelli et al.lll994l) 

The signal from the simulation is also significantly 
smaller than the value predicted using the analytical 
model. As we will see below, the most likely reason for this 
is the fact that the simulation does not include the small 
mass haloes that give most of the signal in the analytical 
case. A limiting factor of the N-body simulations is that by 
construction there is a minimum mass for the haloes. This 
can have a large impact on the predicted average signal as 
smaller haloes are expected to be much more numerous than 
massive ones and they can boost the average signal. In the 
next subsection we explore the range of masses present in 
the simulation. 

4.3 Dependency with the resolution 

In the previous sections, we have shown the results ob- 
tained with the 300 Mpc simulation. In this section we com- 
pare the results obtained with the 50h~^ simulation that has 
a much higher resolution. 

When we compare the mass functions, we find that, as 
expected, the 50 Mpc box contains less massive haloes, 
but it has many more small haloes. A halo must contain 
of the order of 20 particles to be considered a halo. There- 
fore the minimum mass of a halo depends on the resolution 
of the simulation. On the other hand, the maximum mass 
of a halo depends on the volume of the simulation. Large 
haloes are truncated by the boundary conditions of the sim- 
ulation that suppress the power on scales larger than the 
box side. In other words, there is a minimum fe-mode in the 
Fourier modes which is directly related with the dimension 
of the box. 

We compare the average free-free effect in the 300 h^^ 
Mpc and 50 Mpc boxes. Since the 50 Mpc box is 
6^ times smaller in volume than the 300 Mpc one, we 
renormalize the average free-free to the same volume. As 
expected, due to the presence of smaller haloes in the 50 



5 FREE-FREE FROM A SINGLE MASSIVE 
HALO. A NEW WINDOW FOR CLUSTER 
SCIENCE ? 

In the previous section we have shown how the average 
contribution of the massive haloes (groups and clusters) to 
the mean free-free signal is significantly smaller than the 
contribution from the smaller but more numerous low mass 
haloes. In this section we explore the signal of an individual 
halo comparing the prediction from the /3-model with the 
result obtained from the numerical simulation. 

Using the high resolution simulation (50 h^^ Mpc box, 
Afgas = 0.8 X 10^/i~^Mq per particle) presented in Sect. 4, 
we extract the most massive cluster from it in order to com- 
pare its free-free flux with an analytical model. The redshift 
of the simulation is z = 1.6 but the same conclusions can be 
extracted at other redshifts. It is however interesting to ex- 
plore the high redshift regime since the free-free signal could 
be potentially useful to detect clusters in their earlier stages 
of formation and before the gas is too hot to be seen through 
X-rays. In this sense, the free-free emission could extend the 
actual X-ray science in clusters to the range of the radio 
waves. Similarly, the same cluster could be seen through the 
SZ effect but its detection will be harder if the cluster is not 
hot enough. On the contrary, a lower temperature in the 
cluster makes the free-free signal stronger. 

The halo boundaries are defined as i?vir or the ra- 
dius where the over-density equals 200 times the average 
density in the box (according to the common assumption 
of the virial radius r2oo, adopted to our simulation). For 
the most massive halo in our simulation, this radius cor- 
responds to i?vir = 560 kiloparsecs (co-moving) and the 
corresponding total mass of the halo within this radius is 
Mhaio = 9 X 10^^/i-^Mq. This mass corresponds to a group 
of galaxies. We add the flux per particle and compute the 
flux and temperature distortion (see Eqs. ([1} and HQ}). In 
Fig. |6l we show the one-dimensional profile. The solid line 
represents the electron density (in cm"'') as a function of 
the radius. In order to compare this profile with a /?-model, 
the values for Uo and Rc of the /3-model are obtained by fit- 
ting the solid line in Fig. [S]with the analytical profile. The 
best fitting /3-model is also shown in Fig. |6] where the core 
radius corresponds to 1/14 of the virial radius (or p = 14 in 
the notation used above). 

We fix the temperature for the /3-model case to the av- 
erage over the halo particles in the simulation. This average 
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Figure 6. Density profile extraeted from the most massive halo 
in the 50 Mpc simulation box at 2; = 1.6. The solid line 

shows the average electron density profile in concentric shells. 
For comparison, a /3-model is shown (dashed line). The model 
corresponds to a core radius ~ 40 kpc and p = 1/14. The 
smaller box shows the temperature distortion produced by this 
halo as a function of the radius and at u = 1 GHz. The maximum 



distortion is AT/Tcmb 



10- 



corresponds to Thaio ~ 10^ K. With all these ingredients, 
the fluxes for this halo can be derived from the simulation 
and the analytical model. For the simulation case we find 
Ss = 3.67 X IQ-^ Jy, while the /3-model predicts a larger 
flux Sa — 2.80 X 10"* Jy (a factor 7.5 larger). In terms of 
AT, the maximum temperature distortion is about a few 
/xK (at 1 GHz) at the center of the cluster, that is, within 
reach of future planed experiments like SKA. According to 
iBurigana et al.l (|2004l ) . SKA could reach a sensitivity limit 
of 40 njy in one hour of integration and with an angular 
resolution of 1 milliarcsecond in the 4-20 GHz band. More 
massive and denser clusters would produce an even stronger 
signal making the study of free-free emission in clusters at 
radio frequencies an interesting and useful way to study the 
intracluster medium. In Fig. [7] we show a map of the free- 
free signal at 1 GHz in an area containing a more massive 
cluster at redshift 2; = 0.15 extracted from the 300 Mpc sim- 
ulation. In this case, the free-free distortion is of the order of 
1 mK in the cluster regions. At higher frequencies, the tem- 
perature distortion decreases as . That is, at 30 GHz, 
the temperature distortion would be of the order of 1 /iK. 



6 DISCUSSION 

Our results show that there is a signiflcant free- free sig- 
nal at all redshifts up to the time of reionization. Our predic- 
tions are based on analytical models and they are compared 
with N-body simulations. Some assumptions made in our 
model need to be improved, like, for instance, the fact that 
all low mass haloes remain ionized at all times. 

Another important improvement is to substitute the /3- 
model (in the analytical calculations) by a more accurate 
descripti on of the gas in ma s sive h aloes. In particular, the 
model of lAscasibar fc Diegol l|2008l ) assumes a steeper and 
non-isothermal profile for the gas distribution that could 
boost the free-free signal. This model is in better agree- 
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Figure 7. Free-free distortion for a massive halo (M = 6.6 X 
10^'^h~^ Mq) at redshift z = 0.15. The greyscale shows the dis- 
tortion in K and at 1 GHz. The field of view is Ri 40'. The total 
flux in this region is Sg = 2.83 X 10"^ Jy. 



ment with high resolution X-ray profiles in galaxy clusters 
jAscasibar fc D iego 2008; Sanderson et al. 200 ^) and with 
the S Z effect (jPiego fc Ascasibaj l2008l : lOiego fc Partridgel 
I2OIOI ) than the /3-model. 

Another issue that need to be addressed in the future is 
the fact that the free-free effect is significant for a wide range 
of halo masses. This fact, combined with the high range 
of redshifts, makes the computation of the free-free from 
simulations a very demanding task from the computational 
point of view. 



6.1 Comparison with earlier results 

It is interesting to compare our results (based on nu- 
merical and analytical analysis) with those found in the lit- 
erature that use only analytical methods ijHaiman fc Loebl 
1 19971 : [Ohil 19991 : ICoorav fc Furlanettdl2004^ ■ The main differ- 
ence between our analysis and previous ones is that we have 
focused on the bet ter understood regime at low er redshifts 
and higher masses. ICoorav fc Furlanettol (|2004D shows how 
the free-free signal has the maximum contribution at red- 
shift z ^ ?>. The free-free signal is i ntegrated from the be- 
ginning of the reionization [z ~ 12 in lOh|[l999l and 2 ~ 13 in 
Cgoray _fc__Furlanctto 20Qi) until present while we consider 
only a redshift range (0 < z < 7) in which the Universe 
can be considered as fully ionized (on large scales). Also, in 
this work we focus more on massive haloes which are the 
ones that can be considered as fully ionized at all times (for 
2 < 7). In earlier works, only low mass haloes where consid- 
ered in the calculations of the free-free signal. The modeling 
of the low mass haloes is more difficult since they are more 
sensitive to non-linear phenomena. In a future work, we will 
extend our analysis to higher redshifts to include the tran- 
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sition between a neutral and a ionized Universe and a more 
careful modeling of the low mass haloes. 

In[o3 (dii) (see also lOh fc Mackl2003l '). a model is pro- 
posed for the ionizing sources. The model includes the pro- 
duction rates of recombination line photons TVrccomb and ion- 
izing photons A'ion- It makes a clear distinction between viri- 
alized (collapsed) structure that undergo a starburst phase 
and a diffuse gas that is constantly being re-ionized. A halo 
mass function is used to compute the number of active haloes 
(or haloes with a starburst, and ionizing UV flux) and the 
duration of the starburst is set to a constant interval of 
to = 10^ years. Our model is, instead, much more simplistic 
and assumes that all haloes are fully ionized. This assump- 
tion certainly fails in the low mass halo regime. 

In lOhl (|l999l ).the emissivity is computed combining 
an expression for the luminosity of the haloes Li, {M, z) and 
the rate formation of ionizing photons. In this model, the 
temperature is fixed to 10* K. 

In our case, we used a combination of a /3-model plus 
the mass function combined with a scaling law for the tem- 
perature in the analytical case. In the N-body simulation, 
no assumptions are made about the gas profile or its tem- 
perature since these value s are extracted dir ectly from the 
simulation. Recent models (|Fardal et al.ll200ll ) show that the 
gas is seldom heated up to the virial temperature in systems 
with T < 10® K. Instead they are accreted in cold flows (with 
T ~ 10* K). The cold flow mechanism is not implemented in 
our N-body simulations resulting in smaller free-free signal 
from the smallest haloes. 

The N-body simulation includes the contribution from 
both, compact haloes and d iffuse IGM. In the work bv lOhl 
l|l999l ). (and later bv lOh fc Mackil20d3 ). a clear distinction 
is made between the contribution from small ionized haloes 
(that remain ionized for a limited amount of time before 
becoming neutral again) and the diffuse IGM. The authors 
introduce a cutoff flux Sc corresponding to the minimum 
mass able to be ionized and with a temperature of Te = 
10* K. The minimum mass for the ionized haloes with this 
temperature evolves with redshift as Af, ~ 10*(l+z/10)^'^^^ 
/i-^Mq. In our case, the temperature is derived from the 
simulati on and is, in ge neral, larger than the temperature 
used in[oE illii) and ICoorav fc Furlanett"ol (|2004l ). As a 
consequence, our higher temperatures will predict a lower 
free-free signal from haloes. As we mentioned earlier, in a 
future work we plan to include mechanisms such as cooling 
that would reduce the temperature of the haloes (and hence 
would boost the free- free signal). 

The di stor ti on ov er the CMB temperature from haloes 
derived by [o3 (|l999l ) is ATff = 3.4 x 10"=^ K at 2 GHz. 
This result was obtained with no cutoff in the flux of 
point sources (Sc = 0). On the other hand, an estimation 
of the flux coming from the diffuse IGM renders a much 
smaller temperature distortion (ATg = 6.0 x 10~® K at 
2 GH z ) a result later confi r med by ICoorav fc Furlanettd 
l|2004 ) . ICoorav fc Furlanettd |200i) obtain a value of ATg « 
5.0 X 10~^ K for the halo contribution also at 2 GHz. 

Comparing these numbers with our analytical predic- 
tions (see Fig. [2]), we obtain a lower signal at 2 GHz when 
the temperature of the haloes is computed with the scaling 
law Eq. Jg)) (^ 1.73 x 10"^ K). Fixing the temperature to 
~ 10* K the results agree better (r; 7 x 10"^ K). 



7 CONCLUSIONS 

As can be seen in Fig. [S] the mean free-free signal is 
larger for lower mass haloes. On the other hand, large mass 
haloes have larger individual free-free fluxes but they are 
much less abundant. Consequently their contribution to the 
mean signal decreases quickly with increasing mass. 

An interesting result was shown in Fig. [4] where the 
explicit dependency of the average temperature distortion 
with the mass is shown for different redshifts intervals. From 
this plot, it is clear that the average signal is most sensitive 
to halo masses smaller than 10^^ solar masses. Also, from the 
same figure we can conclude that haloes contribute to the 
average free- free signal at all redshifts up to the reionization 
time. 

Even though groups and clusters are expected to con- 
tribute less than less massive haloes to the average signal, it 
should be possible to detect clusters through their free-free 
signal on a one by one basis opening the door for interesting 
studies of the intra-cluster med i um at radio wavelengths. In 
this line, ICoorav fc Furlanettd (|2004l ') discussed the possi- 
bility of de tecting the sign a l from clumps of IG M wi t h AR - 
CADE. In iBurigana etal] l|2004h the model by [o3 (| 19991 ) 
for unresolved free-free emitters has been exploited arriving 
to the indication that the SKA project will be able to detect 
them with deep exposures. 

Future experiments might focus on the detectability of 
individual groups and/or clusters through their free- free sig- 
nature. This signal can be combined with others (SZ, X- 
rays) in multiwavelength studies. 
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APPENDIX A: COMPUTING THE MASS 
FUNCTION 

In ST, the mass function is given in terms of the factor 
V = [(5c/(Ti\/]^, where 5c is the overdensity contrast required 
for the perturbation to collapse, and gm is the rms fluctua- 
tion in the mass scale M. The function 

M^!^^ (Al) 

behaves like an alm ost universal funct ion with respect to 
changes in redshift (|Tinker et al.|[200^ . The quantity p = 
2.775 X lO^^nM/i^M© Mpc~^ is today's average matter den- 
sity. 

This mass function accounts for the fact that the gravi- 
tational collapse of a halo is not exactly spherical but rather 
if follows a triaxial model. For a given cosmological model, 
the evolution of an ellipsoidal perturbation is determined by 
three parameters, namely the eigenvalues of the deformation 
tensor. These are the ellipticity e, the prolateness p and the 
density parameter 5. In their model, the collapse is traced 
independently over the three orthogonal axes and the viri- 
alization of the halo is defined as the time when it collapses 
along the three axes. Since each axis collapses independently 
from the others, collapse along each axis is frozen once it has 
shrunk by some critical factor. 

The term vjiy) in Eq. (|A1|I is parametrized in the ST 
formalism as: 

where v' = av, a — 0.707, and p = 0.3. In the standard 
Press and Schechter mass function (jPress & SchechteJl974l ) 
p — 0. A ^ 0.3220 is the normalization factor given by the 
constraint that the integral of /(;^) in the whole f range must 
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be equal to 1. For comparison, in the original formalism of 
Press and Schechter, the normalization is 1/2. 

In Eq. (|A2|l . the over-density parameter Sc can be 
estimated, given a cosmological model, using the lin- 
ear growth parameter D (z) and 5c{z) (see for example 
iMathiesen fc Evrard|[l99i ): 



5c = 



Do) 
D{z) 



(A3) 



Similarly one could have considered (5c (z) = const and 
o"m(z) — gmD[z)/Do with the same results. The mass 
function, n{M,z), can be easily derived from Eq. (|A1[) and 
Eq. (K2^ . 



We assume a flat ACDM model (Qm + = 1). In 
this case 5c{z) = 1.6866[1 + 0.01256 logm ^miz)]. The linear 
growth factor is given by (|Peebleslll980l ) 



Vx3 + 2 



'3/2/ /3 , r,\-3/2j / 
X ' (x +2) ' dx , 







where x = a/[(l - Q.A)/i2QA)]^^'^ ■ In Eq. 
RMS fluctuation on the mass scale M: 



(A4) 



A2ll (TM is the 



2 1 

'^^^^ 2^ 



dkk'^P{k)W^{kR). 



(A5) 



The window function W{kR) is a top hat function in 
Fourier space 

3(sin(fci?) - kRcos{kR)) 



W{kR) 



(fci?)3 



with R deflned by M = |p_R^. 
can be parametrized as 

P{k) = Ak"T^{k), 



(A6) 

The power spectrum P{k) 



(A7) 



where A is the amplitude and T{k) the transfer func- 
tion. This choice for the amplitude makes it possible to in- 
troduce it in Eq. ([ASP with R = 8h~^ Mpc to obtain the 
value of (jg = 0.8, while the index for primordial power spec- 
trum is set to n = 1 (both values are fro m the fifth year of 
WMAP analvsis. IPunklev. J. et~al] 120091). For the transfer 
function we use the expression given in Bardeen et al.l (|l986l ) 



T(fc) 



ln(l + 2.34i;) 
2.349 



[l-f 3.89g-f (16.1g)^ + (5.46g)^ (A8) 
+(6.71#1-', 



where q = fc(^~^Mpc)/r, and F is the shape parameter of 
the power spectrum (F ~ QmIi). We have compared the 
mass functions obtained with this transfer fun ction and the 
more elaborated one in l|Eisenstein fc Hul l 19981 find ing very 
small differences. For simplicity we use the .Bardeen et al.l 
l|l986l ) transfer function. 



